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In this paper we explore the behavior of the quasi-particle interference pattern (QPI) of scanning tunneling 
microscopy as a function of temperature, T. After insuring a minimal consistency with photoemission, we 
find that the QPI pattern is profoundly sensitive to quasi-particle coherence and that it manifests two energy 
gap scales. The nearly dispersionless QPI pattern above T c is consistent with data on moderately underdoped 
cuprates. To illustrate the important two energy scale physics we present predictions of the QPI-inferred energy 
gaps as a function of T for future experiments on moderately underdoped cuprates. 
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Recently, attention in the field of high temperature super- 
conductivity has turned to characterizing the superconduct- 
ing phase in the underdoped regime. This phase necessar- 
ily differs from that of a conventional d-wave superconduc- 
tor because components of the gap smoothly evolve through 
T c . This leads to a normal state gap or pseudogap above T c . 
Owing to this fact, and unlike a BCS superconductor where 
the order parameter and the excitation gap are identical, there 
are very few ways to probe directly something as fundamen- 
tal as the superconducting order parameter. Within the mod- 
erately underdoped samples, which we consider throughout 
this paper, recent angle resolved photoemission spectroscopy 
(ARPES) experiments have reported JHQl novel signatures of 
superconducting order. Fermi arcs around the d-wave nodes 
above T c rapidly collapse [2] at the transition to form point 
nodes. It was argued that, within the superconducting phase, 
the temperature dependence of the nodal gap has finally pro- 
vided [1] "a direct and unambiguous observation of a single 
particle gap tied to the superconducting transition". A com- 
plementary and equally valuable probe is scanning tunneling 
microscopy (STM) and the related quasi-particle interference 
(QPI) spectroscopy flHSH]. While this probe, like ARPES, 
is generally not phase sensitive, a controversy has arisen as to 
whether these techniques can, as argued experimentally [4, 6], 
or cannot, as summarized theoretically Q5|, LZD > distinguish co- 
herent superconducting order from pseudogap behavior. 

The objective of this Letter is to provide some resolution to 
this controversy by studying the temperature evolution of the 
QPI pattern observed in STM experiments from the supercon- 
ducting ground state into the pseudogap phase at T > T c . We 
employ a microscopically derived preformed pair theory [8] 
that accounts Js}] for all of the complex momentum and tem- 
perature dependence of the ARPES spectral gap as described 
above. The QPI pattern is obtained using the Fourier trans- 
form of the local density of states (LDOS) associated with a 
single impurity. A central result of our study is that, once com- 
patibility with ARPES experiments is incorporated, the ob- 
servation of the so-called "octet model" QPI IIIOL IIHl (called 
Bogoliubov QPI or B-QPI) is a direct signature of coherent 




Figure 1 : Full Green's function in (a) the clean case (thick solid line), 
and (b) in the presence of an impurity (double thin lines) to lowest 
order in Uo, the impurity strength. A thin line represents the non- 
interacting normal state Green's function Go = (u) + £k)~ , and 
the dotted and wavy lines represent the t-matrices t sc and t VB re- 
spectively. Panels (a) and (b) are the basis for our ARPES and QPI 
analysis 



superconducting order: sets of consistent octet model vectors 
can only be found for energies less than the superconducting 
order parameter energy scale. B-QPI does not persist above 
T c , where the order parameter vanishes. 

Our preformed pair theory |8, 9], which has also been suc- 
cessfully applied to the ultracold atoms, is based on BCS-Bose 
Einstein condensation (BEC) crossover theory, which focuses 
on a stronger than BCS attractive interaction and the associ- 
ated short coherence length. The ground state is the usual BCS 
state with a self consistently determined chemical potential. 
As outlined elsewhere [8], we solve coupled equations for the 
fermionic and pair propagators. In the pseudo-gap region be- 
low T*, there are preformed pairs that become non-condensed 
pair excitations of the condensate below T c , For T < T c , both 
non-condensed and condensed pairs co-exist, but there is a 
gradual conversion from non-condensed to condensed pairs 
as temperature is decreased, leading to the conventional BCS 
ground state as T — > 0. 

The pair propagator is related to the t-matrix, t(Q) = 
t S c(Q) + t pg (Q) (with Q defined as a four-vector), where 
A 2 

t sc (Q) = — 7f^S(Q) and t pg (Q) represent the contribution 
to the i-matrix from condensed and non-condensed pairs, re- 
spectively. The non-condensed pair contribution t pg (Q) is 
obtained from a particle-particle t-matrix ||9t] that includes 
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one dressed and one bare Green's function [12]. The con- 
tribution of the non-condensed pairs to the full self energy 
^( K ) = Eq*(<9) G o(<2 - can be wel1 approx- 

imated in terms of a pseudogap parameter A pg (k) because 
tpg(Q) is strongly peaked at small Q for T slightly above T c 
and T < T c . Thus, E(k, u>) consists II 1311 of two terms 
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The two gap parameters, A sc (k) = A sc <y9k and A pff (k) = 
Ap g c/?k, correspond to the superconducting (sc) order param- 
eter and the finite momentum pair gap (pg). The factor ip^ = 
[cos(k x ) — cos(fcj,)]/2 ensures d-wave symmetry. The effec- 
tive gap which appears in the Bogoliubov quasi- particle dis- 
persion is A(k) = A^k with A(T) = Ja*JT) + A 2 pg (T). 

The damping 7 distinguishes the non-condensed pairs from 
the condensate, as was addressed microscopically in detail in 
earlier work lfl4ll . We define the ordered phase through the 
non-vanishing superfiuid density, p s . Microscopic calcula- 
tions IU5I1 based on Ward identities have established that E pff 
effectively cancels in p s so that p s oc A 2 SC as expected. 

The Green's function is G -1 (k, a;) = uj — £k — E(k,w) 
[see Fig. [TJa)]. Here, £k = —2t{cosk x + cosk y ) — 
4t' cos k x cos k y — /i is the normal state tight binding disper- 
sion with t = 300 meV, t'/t = -0.4, and p(T = 0)/t = 
— 1.083. Here, for definiteness, we take T c = 5.9 meV and 
A(T C ) = 52 meV, with n = 0.85, representing a moderately 
underdoped system llrjll . We present the simplest approxima- 
tion ]9J] of our microscopic theory in order to make our calcu- 
lations more accessible. To reasonable accuracy A(T) 
can be found from the BCS gap equation; we approximate [9] 
A pg (T) = A(T)(T/T C ) 3 / 4 forT < T c and A pg (T) = A(T) 
for T > T c . 

The link between ARPES experiments and STM exper- 
iments both below and above T c is based on the com- 
mon spectral function A(k, to) = — Im.G(k, oj)/tt. These 
ARPES experiments constrain the single adjusted parameter 
7 that, nevertheless, has a well understood microscopic ori- 
gin 1 1411 . We find the presence of a perceptible Fermi arc 
(with length > 10% of the Fermi surface) requires 
that 7(T C )/A(T C ) > 0.2. We have importantly determined 
that the behavior of QPI that we report in this Letter is ro- 
bust for the same regime, j(T c ) / A(T C ) > 0.2. Specifically, 
to account for a Fermi arc of about 20% of the length of the 
Fermi surface we take 7(T C ) « 0.5A pg (T c ). The behavior 
of the transport lifetime (see Ref. II 1311 and references therein) 
suggests that 7 has a characteristic cubic T dependence be- 
low T c and linear above: 7 (T) = ~/(T c )(T/T c ) 3 for T < T c , 
7 (T) = j{T c )(T/T c ) forT > T c . Our results are not partic- 
ularly sensitive to the detailed T-dependence, which, for the 
present purposes, could have been ignored. 

In the presence of a non-magnetic impurity, the diagrams 
contributing to the full Green's function up to first order in 
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Figure 2: (Color online) One dimensional cuts for fixed frequency 
ui — — WmeV in the (a) horizontal and (b) diagonal directions. In- 
tensities of peaks fall off rapidly as Cooper pairs lose coherence. In- 
sets specify temperature dependence of q- vector intensity. All octet 
peaks disappear at and above T c , although there is a finite back- 
ground for qi and qi. 



the impurity strength Uq are shown in Fig. |TJb). For a sin- 
gle point-like impurity, the Fourier transform of the first order 
correction to the LDOS is given by 



<5n(q, u>) = -Im 



d 2 k 



(G(k,«)G(k + q,«) 



-F sc (k, oj)F sc (k + q,u)- F pg (k, uj)F pg (k + q, u))] 
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and K = (k, u>). There are two types of FF terms: the usual 
one which depends on F sc associated with the sc condensate 
ifTHl . and a new term (involving F pg ) represented by the last 
diagram on the right hand side of Fig. \T\b); this reflects the 
contribution of pre-formed (or non-condensed) pairs. Impor- 
tantly, this term (which has not appeared in previous work) is 
required for microscopic consistency. 

We have studied the evolution of the QPI pattern with tem- 
perature and find that the QPI pattern remains relatively un- 
changed up to T = 0.5T C . This low T behavior is similar to 
that observed in previous [ 10] T = calculations. However, 
the peak intensity begins to drop rapidly above this temper- 
ature. At T/T c = 0.9 the characteristic (superconducting) 
QPI pattern becomes smeared out because of the increase in 
the number of non-condensed pairs. The octet peaks vanish 
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Figure 3: (Color online) First (Second) row shows diagonal (horizontal) cuts along q x — q y (q y = 0) at u) = —10 meV for increasing T from 
left to right T/T c — 0.1, 0.5, 0.9 and 1.1. The last row plots the B-QPI gap (thick red lines) at the same temperatures inferred from pairs of 
octet vectors along the horizontal and diagonal directions, compared with the ARPES gap (thin blue lines). When the inconsistency between 
the inferred points on the Fermi surface, k x and k y , exceeds 5% the inversion procedure is terminated. The last figure (3o) on the right shows 
the temperature dependence of the energy where the octet model breaks down (red dots), which falls with temperature, like the order parameter 
(black dashed line). The dotted line in the adjacent figure (3n) represents the simple d-wave gap shape. It is clear that, as in experiment [l|] the 
ARPES gap varies smoothly across T c as in a second order phase transition. 



above T c . We find that the QPI pattern above T c is also dis- 
tinctively different from that of a gapless normal phase [E. U I, 
as expected in the presence of a pseudogap. 

In order to elucidate the temperature dependence of the 
QPI pattern in more detail, we plot in Fig. [2] \Sn(q, uj)\ as 
a function of q along the horizontal and diagonal directions 
for u> = —lOmeV. The insets of Fig. [2 a) and (b) show the 
temperature dependence of the peak intensity for the dispers- 
ing vectors. In agreement with the above discussion, we find 
that the peak intensities decrease significantly for T > 0.5T C . 
At T « T c , the peaks associated with the dispersing branches 
qi, (73, (75, and (77 all vanish (the peaks associated with qi and 
(77 merge into the background and thus should not be identi- 
fied as distinctive individual peaks). These insets show that 
the B-QPI pattern is uniquely associated with superconduct- 
ing coherence; this strong signature of T c in QPI is correlated 
with its counterpart in ARPES |01 01 for moderately under- 
doped systems. 

We turn now to detailed plots of the QPI peak dispersion 
that reinforce these conclusions. Our main finding is that the 
dispersive behavior vanishes in the normal state, as is consis- 
tent with experiment [6]. This is shown in Figs. |3a)-(e) and 
Figs.[3lf)-(j) where |<5n(q, uj) | is plotted as a function of uj and 
momentum along the diagonal direction, q x = q y , and hori- 
zontal direction, q y = 0. We identify several non-dispersive 
features in the B-QPI pattern that reflect the normal state. It is 



not clear to what extent the specifics of the non-dispersive fea- 
tures seen experimentally [6] above T c are present here. In the 
diagonal cut at T = 0.1T C and at low frequencies and momen- 
tum q ss 0, we identify the the octet vector (77 [see Fig. |3J a )]-. 
while the peak at momentum q pa 2.5 corresponds to (73. In the 
horizontal cut, we identify the dispersing qi and (75 -branches, 
which for uj — > merge into a single point, q w 2.55. The 
momenta (71, (73, (75 and 97 satisfy the constraints of the octet 
model for sufficiently low uj. 

In the last row of Fig. 3, we probe the effects of super- 
conducting coherence in another way that is strictly confined 
to the T < T c state. Here we display the gap obtained for a 
given temperature from the B-QPI pattern dispersion shown in 
the first two rows, using the octet model inversion procedure 
(0]. Each row was used to generate two sets of points. For 
definiteness, we take a 5% consistency level as a necessary 
condition for the observation of B-QPI. For a given energy u>, 
the inversion procedure yields points on the Fermi surface, k x 
and k y such that A(k x ,k y ) = uj. The ARPES-derived gap 
B91] is plotted for comparison. For frequencies uj > 0.8A SC , 
we find the two sets of Fermi surface points extracted from 
the octet vectors differ by more than 5%. This leads to a fre- 
quency cut-off in the inversion procedure which is plotted in 
the right-most panel of the bottom row as a function of tem- 
perature, where it tracks the superconducting order parameter 
A sc . Consistency between ARPES and B-QPI is found for all 
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T < T c except in the normal state plot in the last row where 
the Fermi arc is clearly seen in the ARPES derived spectral 
gap. One can see from the panel at T/T c = 0.9 that Fermi 
arc physics is also apparent below T c in the distortion of the 
shape of the QPI-inferred gap. Both gaps extrapolate at the 
antinodes to approximately 50meV « A(T). The segments 
of k space where B-QPI inversion is possible correspond to 
the so-called \A\ "Bogoliubov arcs". 

While our results for T <C T c are in agreement with 
those of earlier studies |Iol HI, our results for temperatures 
above T c differ in several important respects from previous 
workA 01 . In contrast to the pseudogap theory proposed in 
Ref. ||7|], we incorporate the presence of Fermi arcs in ARPES 
above T Q . This may explain why the pseudogap QPI pattern 
of Ref. [7] is qualitatively similar to that of the sc state. On 
the other hand, in Ref. [5] the Fermi arcs were incorporated, 
but under the presumption that the gap parameters entering 
the calculation of the QPI pattern and the spectral gaps ob- 
served in ARPES experiments were equivalent. The resultant 
QPI pattern disagreed with the experimental results and it was 
therefore suggested that the pseudogap phase may be unre- 
lated to precursor superconductivity. In contrast, we argue in 
our preformed pair scenario that, while the same formalism 
must be used to address ARPES and STM experiments (as in 
Fig- [D> the ARPES spectral gap should not be directly em- 
ployed in the calculations of the QPI pattern. 

In summary, in our QPI studies two distinct scales appear, 
A sc (k) and A(k), which have the same k dependence. It ap- 
pears naturally in these STM studies that A SC (T), which van- 
ishes at T c , is more influential in the nodal regime, as is widely 
inferred to be the case fljj, |4f| . By contrast A is roughly T in- 
dependent. We extract A SC (T) from the "k-space extinction" 
of the QPI-inferred gap. (As in Reference 11811 . we do not find 
this k-space- extinction point to be directly related to the mag- 
netic zone boundary, as claimed elsewhere 14J].) We extract 
A from the Bogoliubov quasi-particle dispersion, as obtained 
from the B-QPI inversion procedure. A difference with Ref- 
erence [4] is that we find the inversion procedure works well 
even down to the precise nodal point. The failure of this pro- 
cedure when experimentally implemented near the nodes has 
not been understood. 

An important feature of the present theory is its unique abil- 
ity to address the ordered phase at general T in the presence of 
a pseudogap. Other theories are confined to T w or T > T c . 
Because of this pseudogap it is difficult, but important for the- 
orists to find ways to help identify the very fundamental su- 
perconducting order parameter, as we do here. While our mi- 
croscopically based theory JH contains one adjusted param- 
eter Qy(T c )) our theory should not be viewed as principally 
phenomenological. We have found an important correlation 
which makes our predictions very robust. The condition on 
7(T C ) for the appearance of perceptible Fermi arcs in photoe- 
mission is essentially the same as the condition for the dis- 
appearance of octet peaks in the normal state. In heavily un- 
derdoped cuprates (not considered here) we cannot rule out 



a small contribution of B-QPI in the near vicinity but above 
T c presumably reflecting the fact that in ARPES the Fermi 
arcs are not present [18] until substantially above the nomi- 
nal T c . For these more metallurgically complex cuprates near 
the insulating phase, more detailed experiments and theory are 
needed to probe the correlations between B-QPI and ARPES. 

Note added. After this work was submitted we learned of 
related STM experiments but on heavily underdoped cuprates 
by Pushp era/, ill . 
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